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Distant field BHB stars and the mass of the Galaxy II: 
Photometry and spectroscopy of UKST candidates 
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ABSTRACT 

This is the second in a series of papers presenting a new calculation of the mass of the 
Galaxy based on radial velocities and distances for a sample of faint 16 < B < 21.3 
field blue horizontal-branch (BHB) stars. We present accurate BV CCD photometry 
and spectra for 142 candidate A-type stars selected from ubjr photometry of UK 
Schmidt telescope plates in six high-Galactic-latitude fields. Classification of these 
candidates produces a sample of 60 BHB stars at distances of 11-52 kpc from the Sun 
(mean 28 kpc), with heliocentric line-of-sight velocities accurate to 15 kms -1 , and 
distance errors < 10%. We provide a summary table listing coordinates and velocities 
of these stars. The measured dispersion of the radial component of the Galactocentric 
velocity for this sample is 108 ± 10 km s _1 , in agreement with a recent study of the 
distant halo by Sirko and coworkers. Measurements of the Ca II K line indicate that 
nearly all the stars are metal-poor with a mean [Fe/H] = -1.8 with dispersion 0.5. 
Subsequent papers will describe a second survey of BHBs to heliocentric distances 
70 < R < 125 kpc and present a new estimate of the mass of the Galaxy. 
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1 INTRODUCTION 

In contrast with the detailed knowledge we have of the stellar 
content of the halo of the Galaxy, the dark matter content is 
far less well characterised, with both the total mass and the 
size of the halo poorly determined quantities. The accurate 
measurement of the mass profile would provide important 
clues to the nature of the dark matter. For instance, an 
accurate estimate of the mass within a radius of 50 kpc 
is required to establish the fraction of the mass in compact 
objects identified from microlensing experiments (e.g Alcock 
et al. 2000). Quantifying the distribution of dark matter in 
the Galaxy is essential for understanding the assembly of 
the various baryonic components through comparison with 
simulations. This paper is the second in a series of four, 



* Based on observations obtained at the Jacobus Kapteyn Tele- 
scope, the Isaac Newton Telescope, and the William Hcrschcl 
Telescope, La Palma, the Anglo- Australian Telescope, and the 
ANU 2.3-m telescope, Siding Spring Observatory, Australia. 



presenting a new measurement of the mass profile of the 
Galaxy out to the largest radii, r > 100kpc,^ using surveys 
for remote halo blue horizontal branch (BHB) stars. 

The main shortcoming of previous analyses of the mass 
of the Galaxy is the small size of the sample of objects at 
large radii, used as dynamical tracers. Wilkinson and Evans 
(1999, hereafter WE99) calculate the total mass of the Milky 
Way to be M toi = l.9±l% x 1O 12 M , using the full set of 
27 known satellite galaxies and globular clusters at Galac- 
tocentric radii r > 20 kpc (six of which possess measured 
proper motions). This sample must be nearly complete, so 
a new population of distant halo objects must be found in 
order to increase the number of dynamical tracers. Field 
BHB stars are ideal for this purpose. These A-type stars 
are luminous, My = 0.9 (§5), which ensures that they can 
be detected to large distances, and have a small spread in 



t In this paper we use the coordinate r to denote Galactocentric 
distances and the coordinate R to denote heliocentric distances 



© 0000 RAS 



2 L. Clew ley, et al. 



absolute magnitude, so their distances may be determined 
accurately. 

In a recent paper Sakamoto et al. (2003) added new 
kinematic data to the dataset used by WE99, and obtained 
a considerably more precise estimate of the total mass of 
the Galaxy of M to t = 2.5^1,1 x 1O 12 M . Primarily, the 
new kinematic data comprise radial velocities of 412 BHB 
stars, of which 211 have proper motions, at heliocentric dis- 
tances R < lOkpc, with a median distance of ~ 4.5kpc. 
Sakamoto et al. adopted the WE99 mass model, where 
the total mass enclosed within Galactocentric radius r is 
M(r) = M tot /(l + a 2 /r 2 ) 1/2 . Here a is the scale length, 
for which the best-fit value of 225kpc was obtained. For 
r = 20kpc, i.e. a radius containing all the BHB stars used 
in the analysis, this gives M(r) — 0.09M to t- This indicates 
that most of the improved precision derives from kinematic 
data within a radius containing only one tenth of the to- 
tal mass, and therefore relies on extrapolation of a model 
that is tightly constrained only at small radii. In a review 
of mass estimates of the Milky Way halo, Zaritsky (1999) 
emphasises the dangers of such extrapolation (see also Bel- 
lazzini, 2004, for a useful discussion). This motivates a new 
survey for remote BHB stars, at distances that are a sub- 
stantial fraction of the halo scale length a, in order to obtain 
a direct measure of the enclosed mass at large radii. 

Although BHB stars are abundant in the Galaxy halo, 
selection of a clean sample is not straightforward. Samples of 
field A-type stars in the halo are easily identifiable in UBR 
(or equivalent e.g. ugr) multicolour datasets (e.g. Yanny et 
al. 2000). Samples selected by broadband colour include not 
only BHB stars but also stars of main sequence surface grav- 
ity that are some 2 magnitudes less luminous, predominantly 
field blue stragglers (hereafter A/BS), as well as a small pro- 
portion of quasars. The reason that no large sample of re- 
mote r > 30kpc BHB stars has yet been compiled, is that 
the methods developed for separating the two populations 
of A stars (e.g. Kinman, Suntzeff, and Kraft, 1994) require 
signal-to-noise (S/N) ratios that are unfeasibly high for such 
faint stars (B > 18). In the first paper in this series (Clew- 
ley et al. 2002, hereafter Paper I), we developed two meth- 
ods that overcome the difficulties. The methods require rela- 
tively modest telescope resources, yet produce samples with 
high completeness and low contamination. The methods are 
applicable specifically to the classification of halo stars with 
strong Balmer lines, defined by EW > 13A, i.e. A stars 
in the approximate colour range 0.0 < (B — V)o < 0.2. 

Both methods employ a Sersic function fit to the H7 
and H<5 absorption lines. The first method, the Do.is-Colour 
method, plots the average of the width of the two Balmer 
lines against (B — V)o colour. The A/BS stars, having higher 
surface gravity, separate from the BHB stars because of their 
broader Balmer lines. The EW of the CaK line is used to 
filter out a small number of interlopers in the BHB sam- 
ple. We used Monte Carlo methods to establish that with 
(B — V)o colours accurate to 0.03 magnitudes, and spectra 
of S/N — 15A , samples of BHB stars selected by this 
method would be about 87% complete, with a contamina- 
tion of 7% by A/BS stars. Contamination at this level can 
safely be accounted for in the dynamical analysis. Spectra 
of this S/N are also suitable for measurement of the radial 
velocities. The second method, called the Scale width-Shape 
method, plots two parameters of the Sersic fit, the scale 



width b, against the power-law exponent c. The method is 
almost as efficient as the Do.is-Colour method, with 82% 
completeness and 12% contamination, for the same spectro- 
scopic S/N of 15A . Again, the EW of the CaK line is used 
to filter out a small number of interlopers. The advantage of 
the Scale width-Shape method is that colours are not needed. 
For samples of stars with existing accurate {B — V)o colours, 
the first method is preferred, as the contamination is lower, 
and the completeness higher. The gain is small, however, 
and accurate photometry is time consuming. Where accu- 
rate colours are not already available, if telescope resources 
are limited, the best practical solution will be simply to ob- 
tain spectra, and use the Scale width-Shape method. 

In this second paper we describe the compilation of a 
sample of BHB stars in the magnitude range 16 < B < 19.5, 
corresponding to heliocentric distances 11 < R < 52 kpc, 
using the two classification methods described above. This 
first survey (hereafter the UKST survey) starts from pho- 
tographic plates taken with the United Kingdom Schmidt 
Telescope (UKST). In the third paper we describe a second 
survey, for fainter BHB stars, 20 < g* < 21.1, corresponding 
to heliocentric distances 70 < R < 125 kpc. This second sur- 
vey starts from the Sloan Digital Sky Survey (SDSS) Early 
Data Release (EDR) (Stoughton et al. 2002). The fourth 
paper will present the dynamical analysis of these two data 
sets, combined with the 27 satellites and globular clusters 
analysed by WE99. 

The layout of the remainder of the paper is as follows. 
In section 2 we describe the selection of the UKST fields, 
the reduction of the plate scans, and the production of the 
multicolour photometric datasets. We provide catalogues of 
a total of 461 objects with the colours of A-type stars, which 
are candidate halo BHB stars. Section 3 covers CCD B — V 
photometry of 280 candidate BHB stars, and section 4 de- 
scribes the spectroscopic observations of 156 candidates. We 
describe the measurement of the metallicity and the radial 
velocity from the spectra. Section 5 reviews the measure- 
ment of distances of BHB stars. In Section 6 we use these 
results and the methods developed in Paper I to classify 
these stars, and we provide a table of distances and radial 
velocities of the 60 confirmed BHB stars. Finally, in Sec- 
tion 7 we provide a summary of the main conclusions of the 
paper. 



2 THE SELECTION OF BHB CANDIDATES 

In this section we explain how the six UKST fields for 
the survey were chosen, and we describe the compilation 
of calibrated three-band photographic photometric datasets 
in each field. We detail the (dereddened) magnitude and 
colour criteria for the selection of candidate BHB stars from 
these datasets, and provide complete lists of the candidates 
for each field. These candidate lists are suitable, for exam- 
ple, for quantifying number counts of A-type stars in dif- 
ferent directions through the Galaxy. On the other hand, 
the sub-set of objects for which CCD photometry and spec- 
troscopy were obtained was not selected randomly from the 
candidate lists. Rather, both magnitude and colour crite- 
ria were applied. For example, for the CCD photometry, 
if the seeing was poor, brighter objects were selected for 
observation, while spectroscopy was limited almost exclu- 
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b 


RA 


Dec 








J 2000 


SGP 


250 


-89 


55 


-27 47 


SA94 


175 


-50 


2 53 


12 


F358 


236 


-54 


3 38 


-34 50 


F854 


244 


45 


10 23 


-0 15 


F789 


299 


58 


12 43 


-5 16 


MT 


37 


-51 


22 06 


-18 39 



Table 1. Galactic and equatorial coordinates of the centres of six 
fields observed in the survey, ordered by Right Ascension. 

sively to objects with confirmed CCD colours in the range 
—0.05 < (B — V)o < 0.25, since the classification meth- 
ods work well only within this range. In fact, because the 
colour selection criteria were not finalised at the outset of 
the survey, a small number of objects that fall outside the 
selection criteria used to define the complete candidate sam- 
ples, were also included in the CCD photometry and spec- 
troscopy campaign. The consequence of this is that the final 
list of confirmed BHB stars is not a random sub-set of the 
complete candidate lists, and so cannot be used, for exam- 
ple, for quantifying the density profile of BHB stars in the 
halo. On the other hand, since the stars have been selected 
without reference to kinematic information, each confirmed 
BHB star randomly samples the radial-velocity distribution 
function at that distance, and so the sample is suitable for 
dynamical analysis. 

2.1 Field selection 

We selected the fields to give good coverage in opposing di- 
rections above and below the Galactic plane, with a wide 
range of Galactic longitudes. Although the kinematic infor- 
mation is limited to radial velocities, because of the offset of 
the Sun from the Galactic centre the widely different lines 
of sight through the Galaxy provide some complementary 
information on the anisotropy of the stellar orbits, and also 
avoid the risk of the whole sample residing in a halo stream, 
such as those associated with the Sagittarius dwarf galaxy 
(e.g. Ibata, Gilmorc & Irwin, 1994, Dohm-Palmer et al. 2001, 
Newberg et al. 2002). The selection of fields was subject 
to the availability of suitable plate material in the UKST 
archive. This restricts the search to negative declinations, 
and to fields containing pairs of high-grade plates in the 
u, bj, and or/r bands. We chose six fields at high Galac- 
tic latitude, \b\ > 45°, to minimise Galactic extinction. The 
coordinates of the centres of the six fields are provided in 
Table 1, and plotted in Fig. 1. 

2.2 Compilation of uncalibrated three-band 
photographic photometric datasets 

Pairs of photographic plates in each of the three bands were 
scanned by the Automated Plate Measuring (APM) machine 
(Kibblewhite et al. 1984) to produce lists of detected ob- 
jects, and their measured properties, including uncalibrated 
fluxes. These lists were then matched and calibrated. 

In four fields we used pairs of plates in the u, bj, and 
or pass-bands. The exceptions were F358 which had a single 
or plate, and SA94 where we used two r plates. The r band 




Figure 1. Galactic coordinate Aitoff projection showing the lo- 
cations of the six survey fields. The fields were chosen to sample 
different directions through the Galaxy, within the constraints im- 
posed by the availability of suitable plate material in the UKST 
archive. High Galactic-latitude fields were chosen to minimise ex- 
tinction. The dotted line plots zero declination. 

is a slightly narrower and redder band than the standard or 
band. The effective area of the scan in each field, approxi- 
mately 25 square degrees, is listed in Table 2. The effective 
area is the actual area scanned, after allowing for regions 
drilled out around bright stars, multiplied by a factor 0.8 to 
account for incompleteness due to plate flaws, satellite trails 
and merging with neighbouring images. The survey covers a 
total effective area of 150.3 square degrees over the six fields. 

In pairing and calibrating the object lists from the 
plates in each field we followed the same reduction steps 
described in detail by Warren at al. (1991) - in fact the 
SGP catalogue described here is the same catalogue used by 
them to search for quasars. In outline, the processing steps 
are as follows. For each detected object the APM measures 
a parameter intensity that scales monotonically and non- 
linearly with flux, as well as several parameters describing 
the light profile and the shape of each source. The variation 
of the light profile as a function of intensity, for stellar ob- 
jects on the plates, is used to calibrate the intensities onto 
a scale that is approximately linear with flux (Bunclark & 
Irwin, 1983). The main purpose of using pairs of plates in 
each band was to reduce the photometric uncertainties by 
averaging the intensities. The response curve of the photo- 
graphic emulsion is not perfectly uniform, but varies across 
each plate (known as field effects) . In matching plates taken 
in the same filter we gridded the area scanned, and cali- 
brated the intensity scale of the second plate to the scale 
of the first plate in each patch i.e. we forced the field ef- 
fects on the second plate to match those on the first plate. 
A similar procedure was adopted in order to match the ob- 
ject catalogues from different passbands. The bj band was 
adopted as the reference passband and any field effects in 
the other photometric bands were forced to match any field 
effects present in the bj band. 

We used the shape and light profile information to clas- 
sify the objects as stellar and non-stellar, combining the in- 
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Field 


Area 
(sq. deg.) 


Reddening 
(E(B - V)) 


Sample bj 
range (red) 


Sample bj 
range (unred) 


No. 
cand. 


SGP 


24.8 


0.018 


18.47 - 19.77 


18.40 - 19.70 


64 


MTF 


25.6 


0.032 


17.63 - 19.63 


17.50 - 19.50 


128 


SA94 


25.7 


0.067 


16.27 - 19.97 


16.00 - 19.70 


72 


F358 


23.6 


0.011 


17.34 - 19.74 


17.30 - 19.70 


28 


F789 


25.3 


0.029 


15.32 - 19.02 


15.20 - 18.90 


114 


F854 


25.3 


0.048 


16.19 - 19.59 


16.00 - 19.40 


55 



Table 2. Data for the six fields in the survey. Listed are the effective area of each field, the mean reddening, (E(B — V)), the range of 
bj magnitudes with and without the reddening correction and the total number of candidate stars in each field. 



formation from all the plates in order to improve the discrim- 
ination. Finally we converted the intensities to a logarithmic 
scale and calibrated to magnitudes using CCD observations 
of several stars on each plate, as described below. 

2.3 Calibration 

For the calibration we made CCD observations in the U , 
B, and R bands of selected fields. We also used the CCD 
BV photometry of the candidates themselves (i.e. all A- 
type stars). We created U and R sequences from the A-star 
BV photometry using the relations: U — B = 0.5(B — V), 
for < (B — V) < 0.2; B - R = 0.02 + 1.51(5 - V) for 
(B — V) < 1.0. The UBR photometry was then converted 
into the u, bj , or, r natural system of the UKST using the 
relations, 

u=U - 0.01 - 0.01 (U — B) for - 1.4 < U - B < 2.0, 
bj- B = 0.01 - 0.18 (B - R) for 0.0 < B - R < 1.6, 
or- R = 0.00 - 0.01 (B - R) for 0.0 < B - R < 1.6, 
r-R = -0.01 - 0.06 (B - R) for 0.0 < B - R < 1.6. 

The number of stars which denned the calibration curve 
in any band was typically 25. The colour equations quoted 
above are either taken from Warren et al. (1991), or they 
were computed using the procedures described there. 

Bearing in mind the earlier discussion of field effects, in 
establishing the photometric errors there are two terms to 
consider. There are random errors, which can be assessed for 
any band from the scatter in a plot of magnitude difference 
between the two plates as a function of magnitude (Warren 
et al. 1991, Fig. 3). There are also systematic errors, due to 
field effects on the reference plate in the bj band, and these 
are the same in all passbands, due to the way in which the 
field effects in all plates were forced to match those on the 
reference bj plate. These systematic errors do not contribute 
to the uncertainties in the colours. They are larger than 
the random errors over the magnitude range of interest and 
were quantified by measuring the scatter in the calibration 
curves. Table 3 summarises the photometric uncertainties in 
the photographic catalogues. 

2.4 Candidate selection 

Figure 2 shows two-colour u — bj, bj — or /r diagrams for the 
stellar objects in each field, and the selection boxes which de- 





a(bj) 


<t(u — bj) 


a(6 3 - or) 


16 


0.11 


0.05 


0.05 


17 


0.11 


0.05 


0.05 


18 


0.11 


0.05 


0.05 


19 


0.11 


0.06 


0.06 


20 


0.11 


0.09 


0.10 



Table 3. Uncertainties for the photographic photometric data as 
a function of bj magnitude. The reason why the magnitude errors 
are larger than the colour errors is explained in the text. 

fine the samples of candidate BHB stars. This box is slightly 
different for each field because we have allowed for Galactic 
reddening, since the targets are at large distances. In consid- 
ering candidate selection, the r band is treated as identical 
to the or band. This is because A stars have colours close to 
zero, so the difference between the r and or magnitudes for 
A stars, i.e. the product of a small colour and a small colour 
term, may be neglected. Similarly the the difference in the 
extinction correction is negligible, because the reddening in 
these fields is small (Table 2), and the wavelength differ- 
ence between the bands is also small. In this subsection all 
equations involving the or band apply also to the r band. 

The selection box is defined by the unreddened colours 
as follows: 

-0.15 < (u-bj)o < 0.25, 
-0.05 < {bj - or) < 0.45, 
(u - bj) - (bj - or) > -0.40. 

The (u—bj)o cut eliminates objects with large ultraviolet ex- 
cesses from the sample, such as quasars or white dwarfs. The 
(bj — or)o cut minimises contamination by F stars. An addi- 
tional diagonal cut is made because the colour box lies clos- 
est to the stellar sequence in this corner. Converting bj — or 
to B - V using B — V = -0.02 + 0.80(6j - or) the corre- 
sponding range in B - V is -0.06 < B - V < (0.18 - 0.34), 
where the range in the red limit is a consequence of the diag- 
onal cut (seen in Figure 2). Since we are interested in BHB 
stars, which have 0.0 < B — V < 0.20, the incompleteness 
introduced by the diagonal cut itself will be negligible. 

For each field we shifted the selection box along the 
reddening vector by an amount appropriate for the value of 
E(B — V). The mean reddening values for each field, (E(B — 
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Figure 2. Two-colour diagrams for the stellar targets for the six fields listed in Table 1. The selection box defines the complete sample 
of candidate BHB stars. The number of stars in each box, their magnitude range and reddening correction are listed in Table 2. 
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Star R.A. (J2000.0) Dec. (J2000.0) f> 3 u - bj bj - or a(u - bj) a(bj - or) 



CSGP-001 


01 05 28.20 


-26 14 00.79 


19.40 


0.11 


0.10 


0.07 


0.08 


CSGP-002 


01 03 35.32 


-25 54 49.64 


19.05 


-0.03 


0.11 


0.06 


0.06 


CSGP-003 


01 02 44.21 


-26 27 46.75 


18.87 


0.16 


0.34 


0.06 


0.06 


CSGP-004 


01 02 35.56 


-26 12 39.58 


18.52 


0.08 


0.45 


0.06 


0.06 


CSGP-005 


01 01 56.71 


-25 38 03.13 


18.63 


0.16 


0.19 


0.06 


0.06 


CSGP-006 


01 01 44.02 


-26 58 32.59 


19.55 


0.07 


0.31 


0.08 


0.08 


CSGP-007 


00 59 23.47 


-25 48 36.96 


19.56 


0.20 


0.22 


0.08 


0.08 



Table 4. Photometric catalogue of A-type stars selected from six UKST fields shown in Figure 2. The photometric uncertainties in these 
estimates are interpolated from those in Table 3. This table is presented in its entirety in the electronic edition of Monthly Notices. A 
portion is shown here for guidance regarding its content and form. 



V)}, are provided in Table 2, measured from the extinction 
maps of Schlegel, Finkbeiner & Davis (1998). We used the 
reddening relations E(u — bj) = 1A0E(B — V), and E(bj — 
or) = 1.36E(B— V) (note that the reddening vector is nearly 
diagonal) . The relation between extinction and reddening for 
the bj band is A(bj) = 4.035£(B — V). The magnitude range 
of the candidates is different for each field and is listed in 
Table 2 together with the number of candidates. The total 
number of candidates is 461. A catalogue of the candidate A- 
type stars containing their bj magnitude, bj — or and u — bj 
colours and their associated errors is provided in Table 4. 
The stars in this sample are prefixed with a 'O, standing 
for 'complete' sample. 

The photometric calibration observations were, in fact, 
obtained as the survey progressed, and the candidate BHB 
stars had to be selected before the calibration stage. This 
was possible because the clump of A-type stars is clearly 
defined, and because the APM intensities are nearly linear 
with flux. The consequence of this is that some stars ob- 
served photometrically or spectroscopically fall just outside 
the selection boxes. 

In F854 we noted a significant faint over-density of 
fainter (bj — 20.2) BHB candidate stars in the direction 
I — 244° , b = 42° . We later realised these mark the location 
of the dwarf spheroidal galaxy Sextans, which resides at a 
heliocentric distance of 83 kpc. In fact the magnitude limit 
of our complete sample in this field is 19.59, so the complete 
sample is not contaminated by Sextans. But it highlights the 
effectiveness of using Schmidt plates to discover BHB stars 
out to at least 80 kpc. 



3 PHOTOMETRY 
3.1 Observations 

The Do.i^-Colour method of classification requires B — V 
colours accurate to 0.03 magnitudes, so the photographic 
photometry of Table 4 is insufficient. To achieve a higher 
degree of photometric accuracy, BV CCD observations were 
made of candidate BHB stars drawn from Table 4. We 
made a total of 421 observations of 280 stars over 47 photo- 
metric nights during the period 1998 to 2001. As a rule, 
the redder objects in the selection boxes shown in Fig. 
2 were given lower priority for observation, since this re- 
gion will contain fewer objects in the desired colour range 
< (B — V)o < 0.2. The observing dates are listed in Ta- 



ble 5. Three telescopes were used for the photometry, (i) 
The 1.0-m Jacobus Kapteyn Telescope (JKT) at La Palma. 
We used a TEK1 1024x1024 CCD detector for the 1998 ob- 
serving run, which was replaced by a SITe2 2048x2048 CCD 
for subsequent runs. Each chip has a scale of 0.33" pixel -1 , 
(ii) The 2.5-m Isaac Newton Telescope (INT) at La Palma, 
using Chip 4 of the Wide Field Camera, which is an EEV42- 
80 4096x2048 CCD with a scale of 0.33" pixel" 1 , (hi) The 
ANU 2.3-m Telescope at Siding Spring observatory (SSO), 
Australia, using a SITe 1024x1024 CCD with a pixel scale 
of 0.59" pixel" 1 . Generally the brighter stars were observed 
at the JKT and the fainter ones at the SSO and INT. Re- 
peat observations were made between telescopes, in order to 
establish the external errors. 

Throughout each night we observed typically 30 stan- 
dard stars from Landolt (1992). The standards were cho- 
sen to span a large range in colour, and were observed over 
a wide range of airmass. Extinction and colour terms were 
measured every night. The target A-type stars were observed 
as near to the meridian as possible, typically at airmasses < 
1.3. The integration times for each exposure were adjusted, 
dependent on the brightness of the star, the brightness of 
the sky, and the seeing, with the aim of achieving uniform 
S/N for all objects. 



3.2 Data reduction, analysis, and results 

Data reduction techniques employed standard routines avail- 
able in IRAF (V2.11)t for bias, flat-field correction, and the 
removal of cosmic rays. The data taken on the ANU 2.3-m 
SSO telescope was reduced and analysed using FIGARO. 
Twilight sky exposures were used to flat-field the data, and 
a correction was then applied to account for the large-scale 
gradient in the sky. The correction frames were created 
by combining and smoothing uncrowded fiat-fielded image 
frames, from which the stars were a— clipped. This proce- 
dure will, in fact, only be correct if any scattered light in the 
instrument illuminates the CCD uniformly. With hindsight, 
it would have been best to quantify the scattered light, us- 



t IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc. under cooperative agreement with 
the National Science Foundation. 
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Telescope 


Date 




N 


Fields 


sso 


1998 


Sep 13-16 


4 


SGP 


JKT 


1998 


Oct 15, 17-20 


5 


SA94, MTF 


SSO 


1999 


Mar 12-18 


7 


F854, F789, SGP 


JKT 


1999 


April 4, 7-9 


4 


F854, F789 


SSO 


1999 


Oct 11-17 


7 


SGP 


JKT 


1999 


Oct 30-31, Sep 1, 3-5 


6 


SA94, MTF 


SSO 


1999 


Nov 17-20 


4 


F358 


SSO 


2000 


Feb 2-4 


4 


F854, F789 


JKT 


2000 


Mar 23, 26, 28 


3 


F854, F789 


INT 


2001 


Mar 29-30, Apr 1 


3 


SA94, MTF 



Range in V magnitude 


N obs 


-^V stars 




a(B - V) 


15 < V < 16 


24 


6 


0.021 


0.025 


16 < V < 17 


42 


20 


0.021 


0.021 


17 < V < 18 


23 


8 


0.025 


0.031 


18 < V < 19 


85 


38 


0.029 


0.037 


19 < V < 20 


41 


20 


0.040 


0.042 



Table 6. The standard deviation in the value of a single obser- 
vation of V and B — V, as a function of V, evaluated from repeat 
observations in the survey. There are 215 repeat observations of 
92 stars. 



Table 5. The dates and total number of photometric nights 'N' 
for the photometric survey in each of the six fields presented in Ta- 
bic 1. Observations were made using the 1.0-m Jacobus Kaptcyn 
Telescope (JKT), the 2.5-m Isaac Newton Telescope (INT) and 
the 2.3-m Telescope at Siding Spring Observatories (SSO). 



ing a procedure such as that suggested by Manfroid, Selman, 
and Jones (2001). Systematic errors are discussed below. 

Stellar aperture photometry was performed using the 
APPHOT package in IRAF. For the standards we chose an 
aperture radius equal to 6.5 times the measured image full- 
width-half- maximum (FWHM). The seeing was typically 
1.2" - 1.6", giving an aperture radius of ~ 7.8" - 10.4". 
The sky level was measured in an annulus between radii of 
10 and 20 FWHM. We solved for the nightly B and V extinc- 
tion and colour coefficients. The typical RMS scatter about 
the best-fit relations in both B and V filters was 0.015 mag- 
nitudes. This is noticeably larger than the Poisson errors, 
and the excess scatter could be due to scattered light. 

For the survey stars we selected an aperture radius of 
2.3 times the average FWHM of stars in each frame. This 
small radius is close to optimum in terms of S/N. To correct 
to the large aperture used for the standards, we averaged 
the aperture correction measured for ~ 5 bright, isolated, 
unsaturated stars in each image frame. The uncertainty on 
the aperture correction is not negligible, being typically 0.01 
magnitudes, and was added in quadrature to the Poisson er- 
rors for the survey stars. The resulting BV magnitudes and 
errors are found by inverting the transformation solutions. 
The total internal photometric errors for the programme 
stars with 16 < V < 20 respectively are in the range 0.015 
- 0.030 magnitudes in V and 0.020 - 0.040 in B - V. To ob- 
tain an estimate of the size of any systematic errors, due, for 
example, to scattered light, we have looked at the scatter in 
the photometry of objects observed on different nights, and 
on different telescopes. 

The external errors for the photometric measurements 
were determined using the approach of Pearson & Hartley 
(1976) (see also e.g. Kinman, Suntzeff, and Kraft, 1994; here- 
after KSK) , which looks at the range of values R (i.e. the dif- 
ference between the largest and smallest values) in N repeat 
observations of the same star. The estimate of the external 
standard deviation, for a single observation, is then given 
by a ext — R/k, where the coefficient k depends on (Table 
22 of Pearson & Hartley, 1976). Assuming all stars in the 
same magnitude range have similar errors, we can average 
the estimate of a. For measurements of n stars the estimate 
of the external standard deviation is given by, 



_ 1 \ - R, m 
- n 2^ k^N~) W 

i=l,n 

where ki(Ni) is the k factor for star i with iV, repeat 
observations. 

We used this approach to compute the errors on V 
and B — V , in 5 magnitude bins. Note that the error on 
B — V may be smaller than the quadrature sum of the er- 
rors on B and V, since the BV observations were made 
quasi-simultaneously, and on closely the same area of the 
chip, meaning that any systematic errors on the B and V 
magnitudes may be correlated. Before computing the errors 
we have to be confident the sample does not contain vari- 
able stars, for example RR Lyraes. Previous investigations 
of large samples of distant A-type stars contain a significant 
number of variable stars, particularly in the instability strip 
at 0.2 < (B - V) < 0.4. For example, KSK discovered 28 
variable stars from a complete sample of 213 A-type stars, 
and Wilhelm et al. 1999a discovered 56 out of a complete 
sample of 1000. These studies found very few variable stars 
in the region 0.0 < (B — V)o < 0.2 (KSK find 3 stars, Wil- 
helm et al. 1999a only 2). Accordingly, to compute the errors 
we limited our analysis to this colour range. Over this colour 
range we have 92 stars with repeat observations, with a total 
of 215 observations. Table 6 provides the computed standard 
deviations, for a single observation, for V and B — V, in the 
5 magnitude bins. The errors in V computed in this way 
are about 30% larger than the random errors, and indicate 
systematic errors contributing to the overall error budget at 
a level comparable to the random errors. On the other hand 
the errors in B — V are only slightly larger than the random 
errors. 

We have adopted the values in Table 6 as our error 
estimates, reducing these values by vN, for N repeat ob- 
servations. Examination of Table 6 reveals that, in order 
to achieve standard errors in B — V of < 0.03 magnitudes, 
at least 2 observations are required for stars fainter than 
V = 18. This was achieved for the majority, although not 
all, of our targets. 

The photometric survey resulted in BV observations of 
280 stars. The results are provided in Table 7, which for each 
star lists the name, the equatorial coordinates, the mean un- 
extinguished V magnitude, and unreddened (B — V)o colour, 
and the uncertainties for these two quantities. Also provided 
are the value of E(£? — V) used for the extinction correc- 
tions, taken from Schlegel et al. (1998), and N the number 
of repeat observations. We discovered 3 stars with colours 
0.2 < (B — V)o < 0.4 that have large residuals (a^ B -v) ~ 
0.1). These are candidate variable stars, and have been ex- 
cluded from the classification process. These three suspected 
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Star 


R. 


A. (J2000.0) 


Dec. (J2000.0) 


V 


(B - V) 


a(V) 


a(B - V) 


E(B - V) 


N 




01 


02 44.2 


-26 27 46 8 


18.69 


0.22 


0.021 


0.026 


0.016 


2 


CSGP-005 


01 


01 56.7 


-25 38 03.1 


18.53 


0.17 


0.017 


0.021 


0.028 


3 


CSGP-006 


01 


01 44.0 


-26 58 32.6 


19.36 


0.17 


0.028 


0.030 


0.018 


2 


CSGP-009 


00 


56 03.2 


-26 27 30.3 


19.13 


0.28 


0.028 


0.030 


0.019 


2 


CSGP-011 


00 


55 40.7 


-26 53 01.8 


18.65 


0.13 


0.021 


0.026 


0.018 


2 


CSGP-017 


00 


47 42.5 


-25 47 17.9 


19.34 


0.12 


0.028 


0.030 


0.015 


2 


CSGP-018 


00 


46 00.8 


-26 19 38.5 


19.36 


0.26 


0.028 


0.030 


0.012 


2 


CSGP-021 


00 


42 41.7 


-26 20 55.6 


19.26 


0.16 


0.028 


0.030 


0.011 


2 


CSGP-029 


00 


58 09.6 


-27 23 15.6 


19.35 


0.11 


0.028 


0.030 


0.020 


2 



Table 7. Photometric data for the blue horizontal branch star candidates, corrected for extinction. Suspected variables are labelled V. 
Stars that were subsequently found to be quasars are labelled 'q'. Stars in the complete sample have names prefixed 'C', otherwise the 
stars are prefixed T. This table is presented in its entirety in the electronic edition of Monthly Notices. A portion is shown here for 
guidance regarding its content and form. 



Field name 


R.A. 


Dec. 


Nphot 


Nspec 




(J2000.0) 






SGP 


55 


-27 47 


22 


24 


SA94 


2 53 


12 


34 


17 


F358 


3 38 


-34 50 


44 


8 


F854 


10 23 


-0 15 


75 


24 


F789 


12 43 


-5 16 


63 


27 


MT 


22 06 


-18 39 


42 


56 



Table 8. The number of photometric and spectroscopic observa- 
tions, Nphot and Nspec made in each field. The survey contains 
photometric measurements for 280 stars. There are 156 stars with 
spectroscopic measurements. 



variables are marked with a V. The number of stars in each 
field with CCD photometry is listed in Table 8. 



4 SPECTROSCOPY 
4.1 Observations 

We obtained medium resolution optical spectra of 156 sur- 
vey stars over a total of 21 clear nights at the 4.2-m William 
Herschel Telescope (WHT), La Palma, and the 3.9-m Anglo- 
Australian Telescope (AAT), Siding Spring Observatory, 
Australia. The number of stars in each field is provided 
in Table 8. The dates of the observations are listed in Ta- 
ble 9. The aim was to restrict spectroscopy as far as pos- 
sible to stars with confirmed CCD colours in the range 
< (B — V)o < 0.2, accurate to better than 0.03 magni- 
tudes, so that we could apply the Do.i5~Colour classifier. 
However, because the photometric and spectroscopic sur- 
veys proceeded simultaneously, this is not true for all the 
stars with spectra. There are some stars without CCD pho- 
tometry at all, but it is still possible to classify them, using 
the Scale width-Shape classifier. 

In all observations the slit width was adjusted to be 
slightly greater than the image FWHM. Such an approach 
produces a close to optimal balance between minimising 
loss of flux from the target and minimising the signal from 
the sky background, while ensuring that the slit is narrow 
enough that the centroid of the light from the star is known 



to be accurately at the centre of the slit. The latter con- 
sideration is important in order to ensure that systematic 
errors in the radial velocity determinations are minimised. 
Observations in the MTF, SA94, F854, and F789 fields were 
obtained with the WHT using the blue arm of the ISIS spec- 
trograph and an EEV12 2K CCD detector. A 600 line mm -1 
grating was used, resulting in a dispersion of 0.44 A pixel -1 . 
Typical FWHM resolutions were measured from unblended 
comparison arc lines, and were found to be 3.5-4 pixels, or 
1.5-1.8 A. The spectral coverage was 3570 A to 5040 A, 
which includes the H5, H7 and H/3 lines, and the Ca II K 
3934 A line. Observations in the SGP, F358 and MTF fields 
were made at the AAT, using the RGO spectrograph with 
a 1024x1024 TEK CCD chip and a 1200 line mm -1 grating, 
resulting in a dispersion of 0.79 A pixel -1 , and a FWHM 
resolution of typically 4 pixels, or ~ 3.2A. The spectral cov- 
erage of these observations was 3705 A to 4500 A, which 
includes the R5, H7, and Ca II K lines, but not the H/3 line. 

To achieve the minimum continuum S/N ratio of 15 
A , to classify the stars (Paper I), entailed total integra- 
tion times for the targets within the range 800 - 6000 sec, 
depending on the brightness of the target and the observing 
conditions. The longer integration times were broken into 
exposures of maximum length 1800 sec, in order to elimi- 
nate cosmic rays. In six cases the integration was halted after 
the first exposure, when the spectrum revealed the target to 
be a quasar. The stars were observed as near culmination 
as possible, with 85% of spectra obtained at airmass < 1.3. 
Wavelength calibrations were obtained either before or af- 
ter each target exposure. This was sufficient, as instrument 
flexure over 1800 sec. was negligible. 

Nightly observations were made of BHB radial veloc- 
ity standard stars selected from the Globular clusters M13, 
M15, M92 and NGC 288 (see Table 11 and references 
therein). We obtained high S/N observations every night, 
S/N = 40 — 120 A , as templates for cross-correlation mea- 
surement of the radial velocities of the targets. We also un- 
dertook an extensive set of measurements at lower S/N, 
similar to that of the targets, S/N = 15 A , as a means of 
quantifying the radial velocity errors. 

We used standard data reduction routines available 
in IRAF for bias and flat field correction, and removal of 
cosmic rays. We then extracted sky-subtracted variance- 
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Telescope 


Date 


Fields 


AAT 


1993 Sep 12,14,16,17 


MTF, SGP 


AAT 


1994 Sep 7-9 


MTF, SGP 


AAT 


2000 Sep 2-4 


SGP, F358, MTF 


WHT 


2000 Aug 23-28 


SA94, MTF 


WHT 


2001 Feb 27-Mar 3 


F854, F789 



Table 9. Dates of spectroscopic observations, and fields observed. 



weighted one-dimensional spectra using a wide aperture. 
The wavelength-calibration arc spectra were extracted us- 
ing the corresponding traced target-star apertures. A poly- 
nomial of degree three§ was fit to the 30 — 45 arc lines iden- 
tified in each spectrum. 

The measurement of the Balmer line widths and shapes 
followed the procedures of Paper I exactly, averaging the 
results for the HS and lines. After dividing by the poly- 
nomial fit to the continuum, we used our custom Sersic- 
profile fitting software to measure the relevant parameters: 
the line- width at 15% depth Do. is, the scale width b, and the 
shape index c. The routine measures deconvolved quantities, 
and requires the resolution of each spectrum as input. Mea- 
sures of the FWHM from unblended comparison arc lines 
are available for all the object spectra. Comparing individ- 
ual measurements to those with a FWHM fixed at the mean 
value gives absolute differences of at most 3% in any pa- 
rameter. This suggests that the Balmer line measurements 
are relatively insensitive to uncertainties in the estimation 
of the resolution. 

The results of the measurements of the Balmer lines 
are provided in Table 10. Column (1) gives the name of the 
star, and columns (2) to (4) list, respectively, the parameters 
Do. 15, b, and c. The errors on the parameters b and c are 
provided in columns ((5) to (7), in the form of A and B, 
the semi-major and semi-minor axes of the error ellipse in 
the b — c plane, and 9 the orientation of the semi-major 
axis, measured anti-clockwise from the fe-axis. Here the error 
corresponds to the 68% confidence interval for each axis in 
isolation (see Paper I for further details). We also fit the 
profiles of the Balmer lines of the radial velocity standards 
(§4.3), and the results are provided in Table 12. 

4.2 Metallicities from Ca II K lines 

The Ca II K line is the strongest measurable metal line 
present in the wavelength range covered by the spectra, and 
the only useful line in moderate resolution blue spectra for 
measuring metallicity. With the exception of the Am and Ap 
stars, the Ca II K line equivalent width EW, in conjunction 
with the (B — V)o colour, provides a moderately accurate 
measure of the metallicity of the star (e.g. Pier 1983, Beers 
et al. 1992, KSK). 

In Paper I we described our procedures for measuring 
the Ca II K line EW, and for estimating the metallicity. We 
follow the procedures of Paper I here, exactly. The method is 
illustrated in Fig. 3, which plots EWca of the survey stars 

§ in IRAF parlance this is a polynomial of order four 
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Figure 3. Ca II K line EW measurements for the survey stars. 
The curves of constant metallicity for [Fe/H] = -1.0, -2.0 and -3.0 
were computed by Wilhelm et al. (1999a) and the line for [Fc/H] 
= 0.0 is taken from Paper I. Stars to the left of the vertical line at 
(B — V)o = 0.05 have indeterminate metallicities (filled triangles). 
Measurements for the remaining stars are shown as filled circles. 

versus (B — V)o- As explained in Paper I, the metallicity 
is measured by interpolating between the plotted curves of 
fixed metallicity. In Paper I we calibrated the intrinsic accu- 
racy of this method, using high S/N spectra (i.e. negligible 
error in EWca), finding a dispersion of 0.3dex in the mea- 
sured metallicities relative to accurate literature values from 
high-resolution spectra. We have confirmed this value by 
measurements of the BHB radial velocity standards (§4.3) 
used in this paper. For the survey stars, this intrinsic scatter 
of 0.3dex is added in quadrature to the uncertainty in the 
metallicity associated with the accuracy of the measured 
EW. Because the curves crowd together for bluer colours, 
we have not attempted to measure the metallicities of stars 
with colours (B - V)o < 0.05. 

The measured metallicities and their random errors 
(only) are listed in column (8) of Table 10. The majority 
of the survey stars are metal poor, with a mean metallicity 
[Fe/H]= —1.7 with dispersion 0.46, and all but two stars 
have measured metallicities in the range —3.0 < [Fe/H] < 
-0.5. The two metal-rich stars are CF789-093 and CMTF- 
058, which both have measured values close to solar. There 
are 13 stars bluer than (B — V)o = 0.05 (shown as filled 
triangles), for which we cannot estimate the metallicity re- 
liably. 



4.3 Radial velocity determinations 

Radial velocities for the stars were measured by cross- 
correlation of the spectra with our high S/N spectra of 
BHB template stars of known heliocentric radial velocity. 
Relevant details of the template stars, taken from the litera- 
ture, are provided in Table 11. Here the quoted metallicities 
are average values for each globular cluster, and are based 
on high-resolution spectra. We used the XCSAO package 
in IRAF, which uses Fourier methods, based on the method- 
ology set out by Tonry & Davis (1979). After continuum 
subtraction, and transformation to a log wavelength scale, 
the spectra are filtered in Fourier space, to remove low fre- 
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Star D . 15 (j5) b( 7 <5) c(yS) ABO [Fc/H] V d e Classif. 

[A] [A] [kms- 1 ] [kpc] 



CSGP 


003 


31.41 


± 


0.96 


7.25 


0.69 


0.39 


0.028 


1 


519 


-1.4 


± 


0.2 


-46.24 


± 


10.2 


14.5 


± 


2.0 


A/BS 


CSGP 


005 


35.85 


± 


1.39 


9.25 


0.78 


0.55 


0.051 


1 


520 


-1.4 


± 


0.2 


124.13 


± 


14.2 


14.6 


± 


1.3 


A/BS 


CSGP 


006 


35.45 


± 


1.05 


8.40 


0.69 


0.41 


0.029 


1 


524 


-1.6 


± 


0.2 


-103.82 


± 


10.6 


20.6 


± 


1.4 


A/BS 


CSGP 


009 


21.91 


± 


1.15 


5.18 


0.67 


0.38 


0.084 


1 


481 


-0.9 


± 


0.3 


-65.80 


± 


13.5 


45.2 


± 


3.3 


BHB 


CSGP 


Oil 


37.86 


± 


1.07 


9.65 


0.79 


0.44 


0.034 


1 


525 


-1.3 


± 


0.2 


-41.81 


± 


13.2 


17.0 


± 


1.6 


A/BS 


CSGP 


017 


35.13 


± 


1.49 


9.87 


0.93 


0.61 


0.068 


1 


520 


-1.9 


± 


0.4 


27.25 


± 


16.1 


53.5 


± 


4.6 


BHB 


CSGP 


018 


23.92 


± 


1.71 


4.97 


0.61 


0.64 


0.049 


1 


509 


-1.3 


± 


0.2 


89.22 


± 


16.0 


18.8 


± 


1.8 


A/BS 


CSGP 


021 


32.99 


± 


1.58 


8.17 


0.75 


0.65 


0.052 


1 


518 


-2.0 


± 


0.3 


-89.94 


± 


15.2 


52.0 


± 


4.4 


BHB 


CSGP 


029 


31.08 


± 


0.77 


7.27 


0.70 


0.31 


0.024 


1 


518 


-1.1 


± 


0.2 


-68.40 


± 


10.1 


50.7 


± 


4.5 


BHB 


CSGP 


030 


34.77 


± 


1.20 


8.87 


0.78 


0.47 


0.043 


1 


522 


-1.8 


± 


0.3 


-45.18 


± 


13.4 


16.8 


± 


1.8 


A/BS 


CSGP 


034 


28.89 


± 


1.20 


7.09 


0.75 


0.45 


0.044 


1 


517 


-1.2 


± 


0.3 


-41.25 


± 


13.6 


51.2 


± 


4.0 


BHB 



Table 10. Spectroscopic data for the horizontal branch star candidates. This table is presented in its entirety in the electronic edition 
of Monthly Notices. A portion is shown here for guidance regarding its content and form. 



quency continuum variations and high frequency noise, and 
then apodized, before cross-correlation. The peak of the 
cross-correlation function was identified in each spectrum 
and fit with a parabola. The position of the peak measures 
the difference in the radial velocities, and the measured ra- 
dial velocity was then corrected to the heliocentric value. 
We checked that using different functional fits to the cross- 
correlation peak produced very similar results. 

The width of the cross-correlation peak provides an es- 
timate of the uncertainty in the radial velocity. These ran- 
dom errors were typically 7 km s _1 for our survey stars. We 
then measured the scatter in repeat measurements of the 
template stars, over different nights, which was 10 km s _1 , 
significantly larger than the random errors for these high 
S/N spectra. This scatter is plausibly due to the star being 
imperfectly centred in the slit, as confirmed by measure- 
ments where we deliberately offset a template star to the 
edge of the slit. Accordingly we have combined this error in 
quadrature with the random errors. The result is a typical 
accuracy of 12 km s _1 for the measured radial velocity of 
a survey star. The scatter in repeat measures of the radial 
velocity of a number of survey stars, as well as template 
stars observed at lower S/N, is consistent with the errors 
computed in this manner. 



5 DISTANCE MEASUREMENTS 

The absolute magnitude of a BHB star Mv (BHB) depends 
on both colour (i.e. temperature) and metallicity. By assum- 
ing all stars in a globular cluster have the same metallicity, 
the colour dependence may be determined (independent of 
knowledge of the distance to the cluster) by plotting the dif- 
ference between the apparent magnitude of the stars along 
the horizontal branch, and the apparent magnitude of the 
RR Lyraes in the cluster. Stacking the results for several 
clusters, Preston et al. (1991) derived a cubic relation in 
(B — V)o, which we have adopted. To determined the abso- 
lute magnitude of a BHB star of given metallicity, we then 
need to know the dependence on metallicity of the absolute 
magnitude of RR Lyrae stars, parameterised as, 




0.7 0.8 0.! 
Power, c 

Figure 4. Classification of 11 BHB radial velocity stars from Ta- 
ble 12 by (a) the Dq^-CoIout method and (b) the Scale width- 
Shape method. The classification boundaries determined in Paper 
I arc shown as a dashed line. For the Doi^-Colour method we 
have imposed an upper limit of Do. 15 = 36A shown as a solid 
line. One star, M92 IV-17, is misclassificd in the Scale width- 
Shape plot and an adjusted form for the classification boundary 
has been adopted for values of c > 0.85 (solid line). 
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Star 


R.A 


(2000) 


Dec. 


V 


B-V 


E(B-V) 


V 


3 




[Fe/H] 


Source 






















[km 


s - 


l ] 








NGC 288-431 





52 


37.9 


-26 


36 


34.2 


15.22 


0.15 


0.012 


-45. 1± 


1.7 


-1.24 


(2), (3) 


(5) 


NGC 288-421 





52 


38.3 


-26 


37 


03.0 


15.26 


0.12 


0.012 


-39.2 


± 


1.7 


-1.24 


(2), (3) 


(5) 


NOG 288-418 


o 


52 


39.2 


-26 


37 


28.0 


15.22 


0.11 


0.012 


-39.0 


± 


2.9 


-1.24 


(2),(3) 


(5) 


M5 11-78 


15 


18 


27.6 


02 


07 


25.8 


15.02 


0.16 


0.039 


+42.2 


± 


1.1 


-1.33 


(1) 


(5) 


M13 111-58 


16 


41 


37.5 


36 


25 


17.2 


15.07 


0.13 


0.016 


-247.4 


± 


2.2 


-1.56 


(1) 


(5) 


M13 111-26 


16 


41 


30.2 


36 


26 


13.2 


14.99 


0.13 


0.017 


-254.1 


± 


2.0 


-1.56 


(1) 


(5) 


M92 IV-17 


17 


17 


00.9 


43 


13 


04.2 


15.50 


0.02 


0.022 


-127.6 


± 


2.0 


-2.29 


(4) 


(5) 


M92 IV-27 


17 


17 


06.8 


43 


12 


23.9 


15.19 


0.17 


0.021 


-115.7 


± 


2.0 


-2.29 


(4) 


(5) 


M92 XII-9 


17 


17 


30.1 


43 


06 


08.3 


15.09 


0.15 


0.021 


-127.7 


± 


2.0 


-2.29 


(4) 


(5) 


M92 XII- 1 


17 


17 


45.1 


43 


05 


59.0 


15.11 


0.19 


0.021 


-127.6 


± 


2.0 


-2.29 


(4) 


(5) 


M15 IV-44 


21 


30 


02.8 


12 


08 


36.7 


15.79 


0.25 


0.103 


-108.5 


± 


2.1 


-2.22 


(1) 


(5) 



Table 11. Published spectroscopic and photometric data for selected globular cluster horizontal branch stars. 

Sources: (1) Peterson, 1983 (M3, M5 and M13); (2) Peterson, 1985 (NGC 288); (3) Olszewski, Harris & Cantcrna 1984 (CM 
photometry for NGC 288); (4) Cohen & McCarthy 1997 (M92) and (5) Harris 1996. 



Star D 15(75) 6(7*) c( 7 <5) A B 6 d Q N 

[A] [A] [kpc] 



NGC 288-431 


32.43 


± 


0.21 


7.4 


± 


0.3 


8.88 


0.89 


0.16 


0.019 


1 


518 


4 


NGC 288-421 


29.51 


± 


0.18 


7.3 


± 


0.2 


7.44 


0.78 


0.16 


0.017 


1 


515 


5 


NGC 288-418 


31.88 


± 


0.15 


7.2 


± 


0.2 


8.47 


0.84 


0.13 


0.015 


1 


517 


5 


M5 11-78 


30.26 


± 


0.25 


6.8 


± 


0.3 


7.91 


0.82 


0.16 


0.018 


1 


513 


3 


M13 111-58 


31.07 


± 


0.12 


7.0 


± 


0.5 


8.82 


0.94 


0.08 


0.011 


1 


514 


3 


M13 111-26 


32.02 


± 


0.12 


6.8 


± 


0.2 


8.62 


0.86 


0.10 


0.011 


1 


517 


5 


M92 IV-17 


28.89 


± 


0.09 


8.1 


± 


0.2 


8.59 


1.02 


0.10 


0.017 


1 


511 


11 


M92 IV-27 


29.65 


± 


0.09 


8.2 


± 


0.2 


7.73 


0.82 


0.10 


0.011 


1 


512 


10 


M92 XII-9 


27.22 


± 


0.09 


7.7 


± 


0.2 


6.87 


0.78 


0.10 


0.011 


1 


509 


10 


M92 XII- 1 


25.40 


± 


0.08 


8.0 


± 


0.2 


6.18 


0.74 


0.09 


0.010 


1 


506 


11 


M15 IV-44 


27.48 


± 


0.10 


10. 


i ± 0.2 


6.98 


0.78 


0.14 


0.015 


1 


511 


16 



Table 12. Spectroscopic results for selected globular cluster horizontal branch stars. 



M V {RR) = a + /3[Fe/H]. (2) 

The exact form of Equation 2 remains controversial. We 
refer the interested reader to a recent review of this sub- 
ject by Cacciari & Clementini 2003. We use recent obser- 
vations by Clementini et al. (2003) that employ ~ 100 RR 
Lyrae stars located in the Large Magellanic Cloud to de- 
rive (3 = 0.214 ± 0.047 for the slope of this relation. We fix 
the zero point by reference to the measurement by Gould 
& Popowski (1998) of the absolute magnitude of RR Lyrae 
stars, Mv(RR) = 0.77 ± 0.13 mag at [Fe/H] = -1.60, using 
statistical parallaxes derived from Hipparcos observations. 
The final expression for the absolute magnitude of BHB 
stars is then 

Mv(BHB) = 1.552 + 0.214[Fe/H] - 4.423(B - V) 

+17.74(5 - V)l - 35.73(5 - V% (3) 

For typical values (B - V) ~ 0.1, [Fc/H]~ -1.7, this gives 
Mv ~ 0.9. Distances and associated errors are then derived 
using the unextinguished apparent magnitudes, V, and the 
corresponding photometric error. We average repeat obser- 
vations using inverse variance weights, then add the error 
from the systematic error on the metallicity. The result pro- 



duces distance errors of 2-6% for the standard stars and 
6-10% for our survey objects. Distances computed in this 
manner are provided in column (10) of Table 10. For the 
small number of stars with (B — V)o < 0.05, for which we 
are unable to measure the metallicity, we adopted the mean 
value. For interest we also computed distances in this way 
for the radial velocity standards, using the accurate litera- 
ture metallicities (Table 11), and the results are provided in 
Table 12. 

Less is known about the absolute magnitudes of blue 
stragglers and so their distances are much less accurate than 
those of the BHB stars. We have adopted the following rela- 
tion derived by KSK from data for blue stragglers in globular 
clusters published by Sarajedini (1993): 

Mv(BS) = 1.32 + 4.05(B - V) - 0.45[Fe/H]. (4) 

Nevertheless, recent work by Preston & Sneden (2000) 
suggests that field blue straggler stars may have entirely 
different formation mechanisms from those observed in clus- 
ters, owing their origin to mass transfer (McCrea 1964). The 
transfer of mass in a close binary system created by Roche- 
lobe overflow during red giant evolution can produce a sub- 
stantial increase in the main-sequence lifetime of the star. It 
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is certainly possible that we are dealing with a mixture of dif- 
ferent types of stars, with each population characterised by 
distinct absolute magnitude ranges and hence distances. We 
note also that the Preston & Sneden (2000) stars are on aver- 
age redder than the stars dicussed in this paper and may not 
be completely identical. Notwithstanding these caveats, wc 
used the KSK-derived Mv to calculate distances for these 
objects, and stress that the quoted uncertainties reflect only 
measurement errors. We note, however, that Yanny et al. 
(2000) found absolute magnitudes of blue stragglers in the 
halo that were some 2 magnitudes fainter than the BHB 
stars, entirely consistent with the above estimates. 



6 CLASSIFICATION 

6.1 The BHB standard stars 

Before turning our attention to the survey objects, we first 
examine how well our classifications methods work on the 
radial velocity standards observed in M13, M15, M92 and 
NGC 288, which are all known BHB stars. Table 12 lists the 
relevant parameters measured from the spectra. These are 
Do. 15, b, c, and their errors. 

Figure 4 (a) plots D .i5(7<5) (average for the H7 and H5 
lines) versus (B — V)o for the radial velocity standards. The 
dashed line marks the classification boundary determined in 
Paper I, which we here revise slightly to the solid line, by im- 
posing an upper limit Do. 15 = 36A, as explained below. The 
dashed line is denned by the equation Do. 15 = 7918.2(5 — 
l / )(S-2409.4(D-V)g-549.5(B~y)g + 125.2(B-I/) +33.0. 
The plot shows that all 11 radial velocity standards are cor- 
rectly classified as BHB stars (using either curve). 

Brown et al. (2003) have applied our classification meth- 
ods to a sample of halo A-type stars. Their Fig. 6c, sug- 
gests that our classification curve may be a little high, near 
(B — V)o = 0.1 (although note that their plot uses estimated 
rather than measured (B — V)o colours). Referring to Fig. 
5 Paper I, where the classification curve was defined, the 
largest measured value for any BHB star in the KSK sam- 
ple is Do. 15 = 35. lA. Since this parameter is measured to 
an accuracy of ~ lA for our programme stars, by imposing 
an upper limit Do. 15 = 36A, only a very small number of 
BHB stars will be lost from the sample, while contamination 
by blue stragglers (already small) should be significantly re- 
duced. 

Figure 4 (b) plots 6(7<5) versus c(7<5) for the radial ve- 
locity standards. The dashed line marks the classification 
boundary determined in Paper I, which we again here revise 
slightly, to the solid line. In this case, 10 of the 11 stan- 
dard stars lie below the dashed line, and are correctly iden- 
tified as BHB stars. The exception, M92 IV-17, has a value 
of c = 1.02. This is larger than any of the values used in 
the original definition of the classification boundary i.e. the 
dashed classification line is an extrapolation in this region, 
and it appears it may lie too low. Further evidence that our 
classification curve is low in this region comes from the same 
measurements by Brown et al. (2003), referred to above. In 
their Fig. 6d, at large c > 1.0, two confirmed BHB stars 
lie above the curve. Accordingly we have revised the clas- 
sification curve slightly upward in this region. The revised 
classification boundary, shown by the solid curve, is given 
by the expression b = 25.7c 3 - 138.191c 2 + 187.5c - 65.8. 



With this minor revisions we are now in a position to 
classify our programme stars into categories BHB or A/BS. 

6.2 The survey stars 

In Paper I we found halo BHB stars can be separated reli- 
ably from halo blue stragglers provided the spectra have a 
S/N > 15A _1 and EW H7 > 13A (equivalent to the colour 
range < (B — V)o < 0.2). There are 142 such stars in 
this survey. This sample includes 94 stars which also have 
B — V colours accurate to < 0.03 magnitudes. We call this 
sub-sample the 'photometric sample', and use it to compare 
the classifications between the two methods. The remaining 
48 stars contain photometric errors greater than 0.03 magni- 
tudes and so are only classified using the Scale width-Shape 
classification method. 

Figure 5 (upper) plots the 94 stars in the photometric 
sample for the Do. 15 -Colour and Scale width-Shape classifi- 
cation methods. For the Do.i 5 -Colour plot (left) there are 
47 stars below the classification boundary (probable BHB 
stars, plotted solid, hereafter sample A), and 47 stars above 
the classification boundary (probable blue stragglers, plot- 
ted open). In the Scale width-Shape plot (right) there are 
41 stars classified BHB (hereafter sample B) and 53 classi- 
fied A/BS. Given the predicted high completeness and low 
contamination of both methods (Paper I) we would expect 
close agreement between the two methods. Furthermore the 
two classifications are not completely independent, since the 
parameters Do. 15, b, and c are related through the equation 
D0.15 = 2Hn(0.83/0.15) 1/c . 

The agreement between the two methods is indicated 
in the lower plot. On the LHS, the 41 stars of sample B are 
plotted solid. On the RHS, the 47 stars of sample A are plot- 
ted solid. The number below the line plotted solid in each 
of these plots, i.e. the number of stars classified BHB by 
both methods, is 38. Over the range 0.8 < c < 1.05 (lower 
RH plot) the agreement is essentially perfect, with only one 
star in sample B (the single open symbol below the line) 
not in sample A, and all the stars in sample A also in sam- 
ple B (i.e. all the black symbols are below the line). The 
agreement is less good for the region c < 0.8 (corresponding 
approximately to (B — V)o > 0.14, Paper I, Fig. 2), with 
two stars in sample B not in sample A, and with 9 stars in 
sample A not in sample B. This of course is the region which 
makes the dominant contribution to the sample incomplete- 
ness and contamination, and the different classifications may 
be explained by scatter across the selection boundaries. The 
reason for the asymmetry in the numbers scattered (i.e. 2 
and 9) is not clear. 

To reduce the uncertainty in this region, and so min- 
imise the contamination by A/BS stars, we combined the 
classifications in the following way. For each of the 94 stars, 
and for each classifier, we determined the probability the 
star is BHB, p(BHB), by computing the two-dimensional 
pdf in each diagram, given the measured values of the pa- 
rameters and the errors, and measuring the fraction below 
the classification boundary. For each star, we took the av- 
erage value for the two methods, p(BHB), and classified as 
BHB all stars for which p(BHB) > 0.5, i.e. we optimally 
weighted the information in the two plots. With this proce- 
dure, 45 of the 94 stars are classified BHB. Of the three stars 
in sample B not in sample A, 2 are classified BHB. Of the 
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(b-v) 

Power, c 

Figure 5. Classification of the 94 programme stars in the photometric sample (§6.2), using the Do.i^-Colour (left) and the Scale width- 
Shape (right) classification methods. The solid curves are the classification boundaries explained in the text. For the upper plots filled 
circles are stars classified BHB in that plot, i.e. the stars below the classification boundary. In the Do,i$-Colour plot there are 47 stars 
below the boundary and 47 above it. In the Scale width-Shape plot there are 41 stars below the boundary and 53 above it. For the lower 
plots, in each case filled circles mark stars classified BHB by the other classification method. 



9 stars in sample A, not in sample B, 5 are classified BHB. 
Finally of the 48 additional stars only possessing spectra, 
without accurate photometry, 15 are classified BHB by the 
Scale width-Shape method. 

In total there are 60 stars classified as BHB, and Ta- 
ble 13 gives a summary of their kinematic properties. The 
contamination of this sample by A/BS stars will be approx- 
imately 10%. Note that the degree of contamination is not a 
function of distance, since we increased the exposure times 
to maintain an approximately constant level of S/N in our 
sample. Listed in Table 13 are the Galactic coordinates I 
and b, the heliocentric velocity, V©, its la error erv, and 
the heliocentric distance d©, and its la error a&. The last 
two columns provide the Galactocentric radial velocity and 
distance, V ga i and d sa ( respectively. To convert the heliocen- 
tric quantities to Galactocentric quantities the heliocentric 
radial velocities are first corrected for solar motion by as- 
suming a solar peculiar velocity of (U,V,W) = (-9,12,7), 
where U is directed outward from the Galactic Centre, V is 
positive in the direction of Galactic rotation at the position 
of the sun, and W is positive towards the North Galactic 



15 



a 10 

3 



-400 



-200 200 
Velocity (km/s) 



400 



Figure 6. Histograms of heliocentric radial velocities for the BHB 
sample (lines) and the A/BS sample (dots). Bins are 40 km s _1 
in width. 
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Table 13. Summary of kinematic information for the candidates classified as blue horizontal branch stars. 
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Star I b Vq ay cf Q cr d V ga i d ga i 

[°] [°] [km s- 1 ] [km s- 1 ] [kpc] [kpc] [km s" 1 ] [kpc] 



IF358-25 


237.8 


-54.9 


IF854-008 


239.6 


44.4 


IF854-058 


243.9 


41.8 


IMTF-17 


31.6 


-50.6 


IMTF-37 


35.0 


-50.6 



162.4 13.0 43.2 

165.8 12.1 12.1 

173.4 16.8 38.7 

-164.4 16.6 32.1 

-13.3 13.5 42.2 



3.5 40.9 46.3 

1.0 24.6 16.7 

3.7 19.9 42.0 

3.2 -87.8 28.6 

2.9 70.5 38.6 



Table 14. Data summary for the horizontal branch stars contd. 



Pole. We have assumed a circular speed of 220 km s _1 , at 
the Galactocentric radius of the sun (Rq = 8.0 kpc). This 
table is the main result of the paper, and will be used in our 
forthcoming dynamical study. 

Figure 6 shows histograms of the line-of-sight heliocen- 
tric radial velocity of the BHB and A/BS samples. Little dif- 
ference in the spread of heliocentric radial velocities for the 
two samples is evident. A Kolmogorov-Smirnov test gives a 
probability of 0.97 that the blue stragglers and BHB distri- 
butions are drawn from the same population. 

It is of interest to compare the measured velocity dis- 
persion of our sample of 60 BHB stars, mean distance 28 
kpc, with the result of Norris & Hawkins (1991, hereafter 
NH) for a small sample of 9 remote halo BHB stars, mean 
distance 55 kpc. After quadratically subtracting the mea- 
surement errors, in the same manner as NH, the measured 
dispersion of the radial component of the Galactocentric ve- 
locity for our BHB sample is 108 ± 10 km s _1 (the result 
for the blue stragglers is 115 ± 10 km s _1 ). NH measured a 
value of 111 ± 25 km s _1 , in excellent agreement. Our value 
is also consistent with results by Sommer-Larsen et al. 1997 
who find that the velocity dispersion falls from ay = 140 
± 10 km s _1 in the solar neighbourhood, to an asymptotic 
value of 89 ± 19 km s _1 at R > 20 kpc. More recently, Sirko 
et al. 2004a isolate large samples of distant BHB stars using 
the Sloan Digital Sky Survey. They split their sample into a 
bright (g < 18) subsample, which is contaminated by blue 
stragglers at the level of about 10% (i.e. similar to the work 
presented here) and a faint subsample (g > 18), which is 
contaminated at about 25%. If we consider only their bright 
sample then oy = 99.4 ± 4.3 (Sirko ct al. 2004b), at a mean 
distance 16 kpc, again consistent with our value. It is worth 
noting for completeness that the inclusion of their faint sam- 
ple changes their result very little at oy = 101.4 ± 2.8. 



7 SUMMARY 

In this paper we have presented the results of a survey of 
remote halo A- type stars, beginning with catalogues of ubjr 
photometry from UKST plates in two northern and four 
southern high-Galactic-latitude fields. Accurate CCD pho- 
tometry and spectroscopy of candidate A-type stars pro- 
duced a sample of 142 stars with data of suitable quality 
for classification into the classes BHB and A/BS. The final 
sample, Table 14, comprises 60 stars classified BHB, at dis- 
tances of 11 — 52kpc from the Sun (mean distance 28 kpc), 
with heliocentric radial velocities accurate to 15 km s _1 , and 
distance errors < 10%. These stars are suitable for a dynam- 



ical study of the mass distribution in the Milky Way at large 
radii. The stars are all metal poor, with mean [Fe/H]=-1.8 
with dispersion 0.5. The measured dispersion of the radial 
component of the Galactocentric velocity for this sample is 
108 ± 10 km s _1 which is consistent with Sommer-Larsen et 
al. 1997 and Sirko et al. 2004b. The remaining 82 stars are 
classified A/BS. These stars have very similar metallicities 
and velocity dispersion to the BHB stars. 
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